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FQKFJTOSD 

sopite tremendous teohnolo ,ical diseomries and Jovsiopcaante, 
the gas and oil industry of today la stUI plagued with unresolved and 
unoopiained problem and phenomena, especially in the field af maervoir 
engineering. One of these la tha wecroscopie Ilea of fluids in reservoir 
rook fomationa wham gaa and oil am found and produced* Thla thesis la 
a study of bat a snail faoat of tits problaa and has baan conducted with 
tha idea that any work dona In thla im will help con tribute to an in- 
creased knowledge of tha subject* 

tha author wishes to acknowledge tha si mam and helpful sugges- 
tions, advice, and encouragement of Professor Holbrook 0 . Bo toot, Hoad of 
tha Petroleum Knglnoaring 'epsrtaant, University af Pittsburgh, la tha 
projection, pro ;rotaion, and successful acylation of thla study* 
leknowlodgnsat la also aada to thorn authors, investigators, end researchers 
abasa works warn utilised for background and basis of the study without 
which such an undertaking would ham bean Impossible* 
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Petroleum Industry fin; a and obtains its r«s Materials, ,,*s 
and oil, fsr beneath tbs surface of tbs earth In areas of rosk formation 
called reservoirs. These reservoirs are not, as Is usually presumed by 
the uninformed, large voids filled with get and oil which, when tapped, 
say be easily produced. On the contrary, they are formations ef porous 
and pores able reek with email. Macroscopic, open spaces between the 
grains ef the Matrix material in which the gas and oil is stored) s 
factor which makes the production of sueh stored gas and oil difficult 
and, in sons cases, impossible. It is for this reason that the complete 
study of all conditions end probabilities of fluid flew in porous media 
becomes an Important f *ceb In the production of gas and oil. 

tbs purpose ef this study was te investigate and attempt to 
evaluate the inter-relationships ef eons of the majer physical properties 
ef porous reservoir rock which affect the flow of fluids through it. 
Although these inter-relationships have been expressed by several dif- 
ferent formulations, mentioned later, all designed to elucidate fluid 
flew, it was thought that a sera useful end informative modification 
could be developed and applied to the problem of correlating permeability, 
porosity and grain sis# of porous media to fluid flow through such media. 
Such a study, designed to reveal the mechanics ef macroscopic fluid flow 
in porous media, should be useful in evaluating tome of the many factors 
which must ultimately be known if a satisfactory solution of the problea 
is ever to be obtained. Any additional knowledge of these factors should 
contribute to a mere complete understanding of the complexities of gas 
and oil production from the reservoir. 
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to the petroleum engineer, especially on* interested in reservoir 

engineering, a vary important physical property of reservoir rock is 
pe roe ability, which can b* properly defined as a measure of the fluid- 
transit ctluj; capacity of a porous material, or as the ability of a fluid 
to flow within the interconnected pore space* of a porous Material. A 
second physical property, also equally important, ia porosity. In oil and 
gas reservoirs, porosity represents the percentage of total rook volume 
which is available for occupancy by either liquids or gasses or both. In 
connection with the tors porosity, care should be taken to distinguish 
be tween absolute and effective porosity. Absolute porosity is the per- 
centage of void space with respect to total volume. Effective porosity, 
on the other hand, is the percentage of interconnected void specs in the 
total volume. A third physical property of reservoir rock is actual 
grain alas end grain else distribution of the matrix sate rial which forma 
the reservoir rock. 

It has been demonstrated by experiment and enpirieal formulas 
that permeability, effective porosity, and grain else are Interrelated, 
one depending upon the other to varying degrees. In porous material in 
which the pores are inter-connected, there exists no definite relation 
between p* me ability and effective porosity. Ihese, however, are directly 
influenced by grain else, in that grain sire or particle else and its 
attendant factors of uniformity of grain sise, si tape of the jjrain, manner 
of grain packing, and the amount of cenentlft& material between the grains 
(degree of lithi fixation) determine the alee and shape of the pore open- 
ings, the extent of interco neotlons of the pore spaces, and the ratio 
of their volume to the total volume of the rock. 
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Jnve«tdiatione' ,,i ' , ' ,u, ^ J have indicated that actual grain si so 
and effective porosity ass exoallant criteria for Ju. 4 slug permeabilities 
of various porous Media, In amah of tha literature, variations of 
perseabilltj and porosity ad thin individual reservoirs have been found to 
be inter- related to a degree and have been shown by FSttks 1 *' in the case 
of the Bradford Third send of the Bradford field of Pennsylvania to have 
these relationships 

X - 6 ,U9 x 10" l0 ( # 

and K « 5,1 x 1CT 10 | $ ) 6 ’ 9 

where K is the permeability and $ is the porosity. The first formula was 
determined from earlier experiments and the second from later experiments. 
However* these relationships are average and do not apply In some eases or 
throughout the complete possible range of permeability and porosity. Here 
again -rain site, with its attendant factors, directly affects these 
relationships. 

In maklAi determinations of the values of permeability* effective 
porosity and . rain else diameter* the first two are now very easily and 
effectively carried out by ooans of simplified* standardised methods which 
give fairly accurate results. In the case of average grain else diameters, 
this determination bo cornea complicated by the particle itself, in that the 
parti elo Is usually irregular in else and shape, and no knows methods are 
available for defining each a particle in geometric terms. This problem 
beeoaee wore complicated and formidable in the case of an aggregate of Ir- 
regular particles* but has been successfully solved by the use of statisti- 
cal methods for evaluating surface area and volume of the aggregates, 

*^f»iwnc«e are listed in the Bibliography, 







* “X.te-.-l 

**S * **•»* . u • * 






%M * M 



U«l IM]) 



























M M «l |UM 



I « 

4 * 







*1# Vt a to 

Mttth to m air |l 









u 

Mi 





k 



Such sethots ur* baaed un the measurement of tiie so-called average grain 

diameter by relating the irregular parUela to an equivalent geoaetrio 
figure each as a circle or a adhere frora whiefc a dimeter can be calculated* 
Thaee methods determine surface arm and volume free the number and 
weight of the particles, using as a basis the theory of dimensions which 
requires that the surface area be proportional to the square, and the 
volume be proportional to the cube of the dimensions used* 

for the determination of particle also to have physical signifi- 
cance, that is, to relate particle else to definite physical properties 
sush as surface or volume, the arithmetic or geometric mm end the action 
disaster determinations ere superseded by statistical di (tasters. Several 
formulas have been developed and used^ to designate these statleticel 
diameters in tame of various physical properties such as seen disaster 
(length), as an volume disaster, mean surface disaster, moan volu*o-su r face 
dimotsr, and might man disaster* the average grain diameters used in 
this study are man volume diameters, based on the average volume and de- 
fined to be that diameter whose velum divided into the total volume would 
give the number of particles* If d is the average particle diameter end 
f is the total volume, a relationship exists ss 

? e^nd 3 

1 

where n is the number of particles of a given size* If d T *' represents the 
volume ef an average particle, and is the total number of particles 
then 

and 
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II . JttSTORICAL DACK BQUNI) 



Investigation of the flew of fluids tlirouih porous Material at low 
pivasiisM was lid dated by ~«rcy^ In 1356. these studies consisted of a 
series of experiments on the flow of water threw h filter beds free which 
rare/ obtained an empirical formula showing the rate of flow to be very 
nearly proportional to the pressure drop per unit length of porous aadlun. 
In its eloaontaiy for*, this formula, 

v m -(K/}0(flP^ll») 

in ii cates that the flail parae&billty, K, of a porous tutorial proportion- 
ately relates the velocity , v, of flow of a fluid of p viscosity to the 
pressure differential, A P, causing this flow over a length of a 1. Since 
A F is Measured between the outflow and inflow ends, it is negative as in- 
dicated in the foreula above* 

The proportionality constant K is considered to be a specific 
property of tiie porous Material, eepiricaily independent of the dimensions 
of the Material, the pressure differential exerted on the fluid flowing, 
and the viscosity of the fluid. For this reason, K may bo expressed in 
terns of other Measurable physical properties of the porous Material, such 
as porosity and grain sdae. Since the formulation of farcy* ■ equation and 
frost it, his lav, awry investigators and researchers have attempted to show 
relationships of porosity, peroo ability, and grain site* 

In 1360, Seelhel*^ first Introduced grain else into the relation- 
ship by showing that the rate of flow was proportional to the square of 
the average grain diameter. this would indicate a decrease in permeability 
as sand grains decrease in else* 



5 







a & 



•m* iii— m u —i 

a%Sm « w w iw i i i u I in ■■ i ^ *d 







c 



The next investigator to still re train sire In flow relationship 
formal s was Season* who. In IJ?2 «nd 1393, proposed the foraule baaed on 
studies of send filters, 

^ * Rj;( sh/i.) 

where is * permeability constant, h/L is the pressure drop through the 
filter per unit-depth, and dj. is the "effective" dimeter of send particles 
having unifoxwity coefficients of less than 5. "Fffoctivs* grain else is 
defined as the opening which will Just pass ten par cent of the particles 
(by weight) and uniformity coefficient ie defined as the ratio of the sire 
opening which will pass 60 per cent of a sample being screened, to the 
sise which will just pass 10 per cent, then the uniformity coefficient Is 
low, as 5 or leas mentioned above, the particles are acre or less uniform 
in slse, but as the uniformity coefficient increases to higher values, the 
particle sites bocoae widely distributed* 

Is using Karen's equation, If q ie expressed in cubic fret of 
water par day per square foot of filter area, In e&lliaetere, the loss 
ef bead, h, in feat of water, and 1 in feet, then the values of range 
frost 1300 to above U000 depending on the cleanliness of the sen), with 
the usual Halts of 8L for ordinasy sands ranging between 2300 and 3300* 

Jl 

Hasan's equation has been widely accepted and utilised by sanitary engineers 
and Is of considerable historical interest since It represents the first 
attempt to recognise the importance of particle sire and a eethod for its 
satisfactory representation* #avie and Pilsey** retested Hasan's formula 
with the following results * 

q a 2300 (h/l) 

and q « l$0(d^) 2 «^) 6 

uhore T is the temperature of the water in degrees Fahrenheit, is the 
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screen si sc in sdUiiaotsrs through which Just 3U per cent of the particles 
(by weight) pass, and 0 is the porosity intro iuoed into the second fonauls. 

rilichter, w in l$99» published the accounts of his theoretical 
analysis of fluid flow through an ideal homogeneous peeking composed of 
spheres of uni fora else In which he was the first to introduce the effect 
of packing as a factor influencing permeability, which *ave rise to ex- 
pressing the averse pore area in terms of the diameter of the sphere* 

Pitch ter derived a soil fled Poleeuiile relationship in his equation 

q s l0.2Ad 2 h/K^ul 

where q is expressed in cc of water per oeccnd , A is the cross-sectional 
area of the peckin’ In square centime tore, h is the pressure differential 
in cent! tasters of water, d is the disaster of The spheres in cent!** tors, 
i g is a packing constant dependent on porosity, p. is the viscosity of the 
fluid flowing in poises, and 1* is the thickness ulon the direction of 
flow in centimeters. In tome of arcy'e permeability, Glich tor's formula 
would read 

K s 10.2d 2 A a 

For various packings, tabulations have been nade of the different values 
of the peeking constant T, g shown as a function of porosity for porous 
•atari al packed with spheres of equal disaster* 

Kavis and tiisey^ developed a formula 

l/s f a o.oS(#Ao ) 3 * 3 

for use in the determination of the values of Kg where t is the porosity* 

7 

ling recognised the importance of Slichter's formula by usinj it to cal- 
culate the average diameter of irregular particles. It is important to 
note hare that the average diameter of irregular particles determined in 
this manner is the diameter of the particles such that if all of thoa were 
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of this ileoetcr, the packing wool have the same pe raw ability as It 
actually lie* at a glvan temperature and porosity* 

for ths ultimate in percwebiliiy, porosity, and grain alt a ralar- 
tionablps, It remained for ioseny,^ In 1J27, to develop and publish his 
1 sapor Ur. t e»pirical formula 

v « (S^g^F/aS 2 L)(L/l % ) 

where C la a form factor dependent on tha shape ©f tha pores, 0 is porosity, 
g is tha accalo ration constant, A l* is tha pressure differential In grams 
per square centimeter, a is the viscosity of ths fluid flawing in poisss, 

5 Is ths specific surface area in aqua re cent! eaters par cubic cent! :ae ter, 

L is tha thickness of the porous medium in can time tens, 1^ is tha longer 
path of flow in centimeters through tha thickness L, and L/L* is a redac- 
tion factor for tha relationship A P/h. tinea S m S 0 (l - 0) tha above 

equation can be rewritten In terms of S 0 , which is ths specific surface 

»l 

area par unit volume of matrix mater! al, cm , as 

r * <C^ghP/uB 0 *(l * fO*t)(L/L tt ) 

Unaware that Eoaeny had developed tha af oraaen ti one! relationabif. 
Fair and Hatch, ?J in 1933* levs loped an aleoctideRtioal formula by follow- 
lag a Una of reasoning similar to that pursued by Koa any* Ibis aquation 
is 

v - sf Pa h/5u(l - ^) 2 tS o 2 

It can bs saaa that tha factors (C*L/l^)““^ in tha JCoaeny equation hurt 
been replaced by tha constant 5, believed by Fair and Hatch to apply to 
normal, rather compact, unconsolidated porous as t Trial, 

teraan^,!! x tt X9^7$ presented an 5 cp roved derivation of ths 
Koseny and Fair and Hatch aquations* Us in,; as a basis tha general equa- 
tion for fLoe through pipes 
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vnhferw » # 1* the avers, .« velocity, g 1* the acceleration constant, * la the 
hydra xLIc radios, defined as the vaiuaw of fluid la the pipe per surface 
area In contact with the fluid, A P is the pressure differential, kg is a 
shape factor and is the reciprocal of Koaerjy’s font factor C, u is the 
Tiaeoaity of tho fluid flawin >, and is the length of flow path, Carsaan 
set about to substitute or aodify the equation factors so that the equa- 
tion would represent flow of fluid through porous eedia. For v he sub- 
stituted v/jt, later to eodlfy it by the relationship L^/L for How over 
the flew path, JLg, greater titan the L dissension of tte —diusi* Then be 
substituted for s the relationship of #/S 0 (l - ft), and applying, the na'ac- 
lag factor U/lg to the expression A ?A, be arrived at the following 
ressw s a g es— i t 

v * (gjtf 3 A tfc*a - tflHL/L** 

In a series of experiments using a vary wide range of particle 
shapes in unconsolidated nolle. Careen found that the parameter 
of Kosexy’s equation and the parsmter k 0 (L^/i,} 2 of his equation both 
averaged about 5# In the modified Kseewy-Ce«8an equation these parameters 
are expressed as k, with an assigned value of 5, shown by 

r xe a> A P/uS d 2 k(l - tf) 2 l 

12 

fylUe, in his paper on the historical develooo— t of the 
Keseny-Cam— equation, discusses tht ears recent determinations by in- 
vestigators of values greater than 5 for the Koteny-Carsaw constant k 
when the formula has been used in connection with Measurements of flow 
through coos olidated Media* Sines the k shape factor of the Koteny-Caman 
aquation has been held by sost investigators to be constant, the higher 
values of k have been attributed to the (1*A) 2 factor which fuse and 
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Bruce*-* haw called the tortuosity factor, ?. Several deterMnations 
were ®ade by ose end Bruce for the value of k for various consolidated 
media of low perueebility where values were obtained as high as 100* 

Others concur with &oee and Bruce, all indicating that the constant for 
consolidated e»dU varies to higher values and that these higher values 
are attributable to the tortuosity function* 

There has been sons dlsa, xweaent with the concept ef increasing 
values of the Koseny-dartan constant In connection with consolidated media* 
H* C« believes that errors arising from the use of 5 as a value far 

the constant will be far lose than those arising frsm uncertainties in 
other variables, especially the porosity function, expressed as gfi/(l - #} 2 
in the KorenyH arson equation* He has shown, by plotting relative 
resistance versus per cent porosity, that at porosities lower than hO per 
cant there are increasingly different values of relative resistance to 
flow for various equations which all contain a porosity function of aoo» 

id 

sort. allaValle has shown a ainil- r plot of relative permeability versus 
per cent porosity for several different investigators' formulas, all con- 
taining a porosity function* This also indicated that at porosities lower 
than bO per cent the porosity function does net resain constant* It has 
not been possible to test the porosity function /fy(l - 0) 2 In consolidated 
media because it is not readily possible to vary porosity and keep the 
specific surface area constant* 

A recent eoilfl cation has been sad* In the Kesany equation by 
tyllie and Lpanjler* 0 who have combined Site Eoteny equation with properties 
of the capillary pressure desaturation curve to obtain the formula 

K • <y 2 /2*5F 2 f{) J* dS,/P 0 2 

ehare £ is the permeability coefficient, y is the interfacial tension, F 
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la tha formation fester, $ la the porosity fraction, 5 W is the fractional 
arottdng phase saturation, and P Q la tha ca^Haiy pressure, This aquation 
introduces nan physical factors nfdch are eors easily masurebla in 
laboratory deterotn&tions and consequently should further 'tea continuing 
search for aero accurate and useful relationship formulations for fluid 
flow through porous stadia. 
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For the purpose of carrying out this Otway, data of axporluwmtal 
«nd calculated value* for pesmeability, together with aecompanylng valuta 
for porosity and average grain eixe, wn selected front F*ttke’a i#2 
studies of the Bradford Oil Field, Bradford, FWmaylvanla. Thmse data 
had oeen derived firm taat determinations made on representative core 
samples obtained fro* various motions of the field at random depths in 
the Bradford Third sand. The Bradford Third a and is m excellent repre- 
sentative of a consolidated porous medium whose general ehareeterlstlee 
art lorn permeability and porosity and quite fine sraln site with irregular 
chape and else distribution and a moderate to high amount of lithl fleet! ecu 
Sines practically ell oil and gat reservoirs art located in eoasoli dated 
roek formations, s representative sample of this type ef porous medium 
ms ehoesn for application to this study, 

Is s first approach to the problem, an investigation was made ef 

permeability, porosity end grain else relationships then applied to to 

13 

equation modified by the author from one developed by Traxler and Baum 
from the Foiaeuille aquation of flew through pipes. This equation was sir- 
pressed by Traxler and Baum as: 

» f s(32W^^ (1) 

where l> c is the effective pore diameter in centimeters, q is the rate of 
floe ie cubic centime term per second, p. is the viscosity ef the fluid 
flaming in poises, L is the length ef flow path in centimeters, A is the 
cross-sectional stem in square centimeters, ff ie the porosity, and d P it 
the pressure differential in dynes per square centime tor, Jt can be 
readily seen that the function $*L/a£ P of this equation mill be equivalent 
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to tho peree ability coofficiant K, in d&rciee, whan ti» unit* an* changed 
to theee of th* Deny aquation* 7bm $ K m {qlOO r tL/A9.i3692 X 10~? A P) or 
9 *j69& X id " 9 K m q^L/AAf* Substituting in aquation (1) above give*: 

$ 9 m 5.619 X lO^CK/f)* (2) 

Khar* t la the pa r e * a b i li ty coefficient in darelaa. 

fable 1* Appendix I, haa bean presented to shew calculations which 
resulted from tha application of tha pe ratability and porosity data ee lac tad 
for tiiia study to aquation (2) above, aolvinj for I> e , th* effective 
average pore dia.ee tar* Suiting^ haa stated that th* affective average 
pore dimeter la roughly one- fifth of the average grain diem tar* Figure 1, 
Appendix II, show* a plot aad* iron tha data in Ibble I of d, th* average 
grain disaster, versus .9 , the affective average pore diameter, resulting 
in a straight Una ehoee slop*, which represent* the ratio of d/$ 9 , is 
approximately 25# The value of this ratio is considerably greater than 
the one ef 5, (l/5d a %)$ proposed by *uttin> It should be explained 
here that only tha results of th* calculations using Fettles’* 193U data 
were used in plotting Figure 1* the results of tha ealoilatiooa using 
Fettle' s 1933 data produced no conclusive plat* 

Aa a second approach to tha problem, an investigation was made of 
permeability, porosity and grain else relationships whan applied to a 
formulation derived by equating the Koaeny-Caman and the rarsy equations* 
Ibis second study was nad* co tha baeia of tha specific surface area, 

5 0 , ef porous medium par unit volume of matrix sate rial. 8y taking the 
baste Koaeny*Csrmsn equation solved for S 0 £ one obtains i 

s *. — 

0 5(1 - V) 2 <iit 



( 3 ) 
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It cun to aeon that In equation (3) or t ho preceding page the function 
(A A P/q >L) la equivalent to the K of arcy's equation when the unite of 
P and A f are changed from poioeo to centtpolees anl from gra m per 
square centimeter to atsooph- res reepestivsly. Thu*, K « qlOQul/ 

(A9.073 x 10 ^*4?) or 9,673 » 10* 6 £ x (AAP/<*t 1)“ X . Substituting in 



equation (3) elves s 



2 . g l< 3 

° 5(1 - 7 ) 29*673 * iCT^S 



<J»> 



fine* the acceleration constant, g, 1* 930 centimeters per second squared, 
equation (h) abort can new to resolved as 

u 

s 0 X k$0i. jJLy-* (*A) S (5) 

which is the equation used in this second approach to the problem under 
study. 

Table II, A pendix I, has been prepared to show calculations re- 
sulting fro® the application to equation (55 above of the permeability 

20 

and porosity data selected for this study. X t has been shown by Harlan 
that S 0 « Z a /Z v d for Irregular particloa where Z S A V is the ratio of the 
surface to volume factors and d is the average grain disaster. In this 
ratio, Z g x S/d # , where S la the surface area end d^ is the statistical 
particle diawe ter based on surface art a, and Z r « v/d v ^, where ? is the 
volume and d y is the statistical particle diameter based on volume. The 
ratio Z 9 /Zy has been iniicated by Fair and Hateh^ to be a useful measure 
ef particle shape and has been ;ivcn a value of 6.0 for spheres, 6.1 for 
rounded particles, 6.U for worn particles, 7,0 for sharp particles, and 
7.7 for angular particles. Figure 2, Appendix II, shows a plot made of 



the values in Table II of & Q versus l/d which roughly for* straight line* 

with a slope of 143 for line 1 using Fettle 1 * 1934 date and a slope of $? 

for line 2 using Tbttko’s 1933 data, the slope normally should represent 

the values of the ratio ZjZ^ which lie eosewher* between 6.0 and 7.7. 

As s third opps'oach to the problem, an investigation was wads of 

permeability, peeoeity and grain ales relationships when applied to a 

formulation derived by equating the Kosery-Canaan and the arcy equations, 

ss previously wenticnod, except that the fsrmlation was noil fled for the 

use of average grain size in lieu of specific* surface area. Since S Q is 

2 2 

proportional to l/d, S Q can be replaced by l/d in equation (4) above, 
fhore ia also a new constant, k 1 , introduced into the equation which it 
equivalent to the function l/(£ # /Zy) 5 where 5 is the Koxeny-Carasa con- 
stant and ZgA v la the ratio of the eurfaee to voluwe factors. The re- 
salting aquation is expressed ss: 



1 

dT 



*LiJL — 

(1 - 0)*9.675 X KT 6, 



( 6 ) 



uy cubs titu ting 900 centimeters per second squared for g in equation (6) 
end rearranging to solve fbr k*. Use equation can now be expressed ss: 



k» - 




(?) 



Table IZX, Appendix 1, has been prepared fro® calculations derived 
by applying the data used in this study to the factors of equation (7)} 
namely, K/d 2 and 1.013 X loV/(l - 0) 2 , the values of these factors were 
plotted ae K/d 2 versus 1.013 X 10®jtfty(l - 0) 2 > shown by figure 3, 

Appendix II , to be e straight line with s slope of 0.0001052 in the case 
of line 1 using Tettke’s 1934 data and a curve with s variable slope of 
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urJmown value In th« cm* at line 2 using; Fwttka's 1933 data* Line 2 
could not be expressed as a straight Has by rs plotting on seui-lo or 
losing graph paper* tbs slops of 11ns 1 of HpK» 3 la equal to the 
constant k* of aquation (7) on tbs preceding page* 

As previously stated, k* is equivalent to the function l/(Z 8 /s v } 2 5* 
Using the value of ths slope date rained for line 1 of Figure 3 for k*, the 
value of 2 # /2 y as obtained froa the relation ?. # /S T • (lA*5)* is i*3*5* 
Somally, values for 2^/2^ lis between 6*0 and 7*? as previously eontioned. 

In the two preceding investigations, the porosity function, 
rf 3 /(l - flf) 2 , although it varies with the value of porosity, is considered 
a constant factor in that it retains ths same fixed relationship regard* 
lose of the source of data* For a fourth approach to the problem, an in- 
vestigation was aade of the pe inability, porosity and grain alee relation* 
ships when applied to a formulation developed to exprese the porosity func- 
tion aa a variable, dependent upon the source of data, with the peraee* 
bility and groin else being constant factors regardless of the scare# of 
date* Urns, the porosity function would be e constant factor only for a 
particular reservoir formation or even s particular section of ths 
reservoir formation* 

To develop this foraulatlcn, ths basic Koseny-Caraan equation was 
recast to include particle disaster, giving* 

, - k BZ < 3 > 

» l (l - w 

whose units have bean previously defined* Solving for the porosity function 
gives* 
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(1 - 




q '±1 

AAP g k* d 2 



(9) 



Sine* tho function qG/AZi ?, whew expressed In taree of the trey equation* 
it equal to the permeability coefficient K of "arey’a equation, tubeti bu- 
tton of 9*670 X KT 6 .*: for qxL/kH? can be wade in equation (9) above giving 



£ 

d-0} 5 



9*670 I 1<T°S. 
'™fc k» d 2 



( 10 ) 



Assujclng that the constant k* would be reflected In e recasting of the 
porosity function, * new expression of this function can bo substituted 
for k*#fy(l - #) 2 in equation (10) above giving 

(jg) n s l / 1.013 I loV ( 11 ) 

or, by rearrangement, 

I m 1.013 X 10^ d 2 {4) n (12) 

where n Is a constant applicable only to a particular formation or section 
of a fo matt on, or a porous wedluau 

Table IV, appendix 7, has been prepared from calculations derived 
froo the application of Fetttce's 193U data to the functions of equation 
(12) above. Figure 1*, Appendix II, is a leg-leg plot of these functions, 
S/d^A*013 X 10^ versus which fora* a straight line with a slope of 3. 
Since 

n . ^ i/x - ou 1 11)3 
then n a 6, the elope of the straight line. 



ti 




iv. co» mj? Am hbcusshk or ssgulw 



Before discussing the Jesuits, It should bo pointed out her# that 

too value* of po realty data used in this a tody represent absolute 

2 

porosity. However, Fcttke has siwnm that In a majority of eases the 
difference between absolute porosity and a iow*r value for effective 
porosity in the Bradford third sand ia but a natter ef 0*1 per cent 
average. There are a fee isolated cases of a deviation illicitly greater 
than 0*1 per sent up to about 1.0 per cent, for purposes of this study, 
tbs values of effective porosity mere aasueod to be the ease as those for 
absolute porosity* This fact can explain soae of the slight deviations 
of tbs points used to plot the various lines end curves shown by Figures 
1, 2, 3 end u. Other deviations of the points can best be explained by 
the fact that et the tie* the data were generated, the accuracy of the 
determinations ef values for portability, porosity end grain six# was not 
nearly as good as that of ths present day, due to greatly improved and 
standardised eetboda for determining sudd value*. 

the results of the author** first investigation clearly sho* that 
the average pero diameter in a rather fine-grained, consolidated medium 
decreases considerably free the average for unconsolidated media, staled 
by Wotting to be roughly one- ft fth of the average grain disaster* This 
can be best explained by tie non-uni f omi ty of particle else and shape, 
the fine grain size, and by the degree of lithiilcatloa or cotenUng 
found in consolidated sella. It is considered by the author that the 
latter exert* the greatest effect upon average pore disaster since count- 
ing to any extent tenia to restrict, reduce and plug the pore channels of 
the rock* To utilise this esthod of correlation of permeability, porosity 
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and ;r*un sis©, lb would be necesuary to Istersalne the relationship be- 
tween »v?ra;,o effective pore disaster and average grain diameter for each 
particular reservoir or aeotion of the reservoir# 

In the second approach to the probiea, the results of the in- 
vestigation indicate that the value of the bozeny-Caraan constant for 
consolidated media Is ’such higher titan the assigned value of 5* Py way 
of exa spl* f preauao the consolidated aeditw to be coepoeed generally ©f 
•harp particles *here the value of the ratio of surface to volume factors, 
Z 9 /Z v , is 7#0. It has been previously shown that S B * Xg/i’^d# Ueing as 
a selection frox Fettke'a 1934 data, d r .'J093 centime tar, K » #0034 arcy, 
and 0 - .145, it is possible, by resolution of equation (4) of section III 
of this study, to solve for a value of the ieseny-Caraaa constant in lieu 
of the value of 5 noraally used* 4^ is resolved as 7/.0093 or 753# 

Colvin, for the unknown Eseeny-Garaan constant using the above value# in 
Use equation gives a value of approximately 220# by u»in r some of the other 
sets of values for pemsability, porosity, and grain else from the 1934 
data to solve for the Koveny-'laraen constant in the same wanner as described 
above gives values for the constant in the general vicinity ef 22Qtrlth 
suae variations# The results of si ail- r application of Fet‘Jce*e 1933 data 
were not as conclusive, with the values showing a such wider variation, 

30 to 107 , for the iosteny-Caraan constant.# 

As in the second approach, the third approach to the problem also 
indicated that the Soseny-Canean constant for consolidated sails is aunh 
higher than 5# It has been previously stated in Section III of this 
study that the slope of line 1 of Figure 3, Appendix XI , was 0#OOG1C52 and 
equivalent to the function 1 /(£*/£*) k where k is the KeaeryCaraan 
constant. If, for example, it la again assumed that £ a /2 T has a value of 
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7.0 far this mediae, it ia possible to soiv* for * value of other than £ 
for Us* Ka^sny-. arson constant* in this Cioo, the value w-s resolved as 
19k for the con.'itant* Cince a emU, 4# via lion In the drawing of line 1 
of either Figure 2 or 3 result « in a change in the value of tf»e elope of 
the lino, it can be presumed that this fact cause! the dffare. ce in 
values for the £e sen?-"' &rx*« constant dele reined by the tec- Invalidations. 
Therefore, it can also be presumed that in the case of the 19Va data used 
the average value for the ko*eny~C«r®*n constant lie* a ee ewt' . ere between 
19u and 220* Wo stailsr success was achieved in usin.; Ft tike’s 1933 data. 

I t is believed that the 1933 data were no* truly representative of the 
reservoir sand nor were they in thesselvee complete enough to pondt drw 
In* ary conclusions. 

The fourth approach to the preties was an attempt to show the 
porosity function as a variable rather than a constant fsefcsr in lieu of 
a variable for the Koveny-Caruar! constant as shown by the second and third 
Investigations. An an plication of the 193U data to equation (12) of part 
III Indicates that the n of the porosity function, fP, is approximately 
equal to 8, as lets rained fn» the slope of the line plotted in Figure h. 
Appendix II, with higher and lower variations in sosas cases. i incs tim 
porosity function is an arbitrary selection by the author, it follows 
that this particular function, bein? an exponential, could have a variety 
of fores all of which would, when solved numerically, „ive substantially 
the sane values. xf i was adopted because of its alar: 11 city if for no 
ether reason. Although this eetttol shows that the porosity function can 
well be a variable factor, it Is believed that the works of ooe and Pruoo^ 
and others have conclusively shown that the porosity function rc«ai ns 
constant while the koseny-Oaraan constant becomes a variable factor In 






in 



I %* 




m -4 



* •* • 



«9 .ail a* 













• ^aa,U 









*« 

I 441 






19 



21 



con* oil -is. led media. 80 attempt wa» mJ« in the fourth investigation to 
aait* os* of Fsttka’o 1933 data bocs-use of previous Inconsequential 
results. 

la all of the author** four Investigations, it hae bean shown 
that the procedures, equations, ml constants applic bis to tbs satirical 
formulations derived from the studies of unconsolidated porous madia do 
not hold true in ths case of consolidated porous stadia. However, it also 
has bacn shown that for a ^iven or particular consolidated porous rock 
for** lion, the estpirical formulations apply when ths normal constants ars 
changed to those values which fit ths si tcation, as in the case of this 
study of ths ©onaoliJatsd Bradford Third sand* Ths evidence 'oea not 
eonclu 3 ivsly prove whether the porosity function sh-mld be a variable or 
a constant factor in pereseM. iity, porosity and grain size relationship 
formulations. l*uch mo re investigation and research ®ust be carried eat 
alen,; those lines before a conclusion can be reached as to which theory 
is correct, although as mentioned in the precadlng paragraph present 
evidence points to the porosity function being a constant factor* It is 
very possible that both ths porosity function and tha toaeny-Caraan con- 
stant night be found to be variable# in the cose of consolidated parous 



media 



V. C3NCUBXGIB 



The field of *=»I1 parti, els eta tin ties and asasarewenis 1* exten- 
sive, coxples end relatively unknown. For those reasons, the study of 
particle site In relation to portability and porosity of porosis madia is, 

of itself, extensive, complex ani relatively unknown, especially in the 

21 

case of consolidated porous sella, fho uthor o race with Alaska t who 
stated, in effect, that the use of snail pxrticle etatiwtice and weasure- 
aunts in deterednations of permeability was so complex and unresolved that 
ha presently utilized tssthols of deicrs&nin^ permeability and porosity 
far outmatch ary awtfced usitt; ,;r*un size or its attenisnfc properties to 
determine sam, However, the author feels that it is essentially im- 
portant to intensively study the field of particle sis© in order to ;ain 
a thorough knowledge and uniers tending of macroscopic fluid How In porous 
media in the lijh t of ^rain size and its attendant factors of packing, 
uniformity of size and shape, and decree of lithifi cation, all of which 
affect fluid flow in such media. 

Through an Intensive study of this nature, it is possible that 
•one day a feraalation will be developed which will truly hold urrer all 
conditions ani in any circusetancee whereby tbo physical properties of 
permeability, porosity, and grain »i*e with Its attendant factors, can be 
successfully correlated. The four methods of correlation espl-ywd in 
this » tady can be used In a restricted degree when applied to the con- 
solidated porous media of particular reservoir* or of sections of the 
reservoir* only after considerable research to date mine the proper 

c ns tant* applicable In each case. In the saae vein, and of considerable 

22 

interest, is the Method developed by Hyder for saaktn; permeability 
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asaaursaent# without the use of norm samples by mums of an empirically 
i«volopei formula, g'nd only for on* particular sand, which uses only 
istersdrad weights of sieved sttni particles which haw* been passed 
through a series of U. $* Standard s lavas. 

Another possible investigation or wsthod, suggested to, but un- 
tried by the author, would be to plot the result* of sieving data froe 
a cor* aasiple (per cent by weight of sample through and on each screen 
else) versus the average screen sire for that per cent by weight, Free 
the resultant cures, select s s representative" per cent value and deter- 
aln* the corresponding average screen sits or presueed average grain 
disaster. Using equation (2) of Section III of thie study, detoreine 
which per cent by weight value and cormspotvAlrig screen si as or presaeed 
average grain site can be used to satisfy the equation, and then determine 
whether or not this saw* per cent value can be used in connection with 
all core staple data from cores taken froe the s&ae foraatlon or section 
of the foraaUcn. 

A saswuticsn of the results of this study does not reflect the 
original purpose of the study. 1 t w*s thought that an evaluation could 
be aade of the Interrelationships of the physical properties of porous 
ruck, permeability, porosity and grain si so, as mil ss a erfifi cation of 
OJd.atU. 0 fc foreul'tiona used to express these Interrelationships of physical 
properties, hy using data froe highly consolidated por nts walla, it was 
hoped that If a -soil fleet! on wore developed for effect! vsly expressing ths 
Interrelationship of perseability, porosity, and grain site, it would hold 
unlar any cl ream tances and conditions for any type aeilus, either con- 
solidated or unconsolidated. However, the complexities of grain site 
study, the ieerth of sufficient and accurate data, and the apparent non- 
conformity of application of aoii fled f emulations to consolidated porous 
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media data made such enleavort impractical enl Impossible* 

It le believud that the results of this study, although meager, 
can contribute some benefits in a better under* tending of the problem* 
the ratio of pore voluac to grain else liaaoter has been shown to decrease 
markedly with an increase in cementation and a corresponding decrease in 
grain else such ae found in a consolidated porous media like the Bradford 
Third sand* The value of the Sose«xr»Carsum constant has been shown to 
increase to very high values in consolidated porous media* This is prob- 
ably due to the large increase in values of the tortuosity factor which 
becomes greater as the rook becomes aore consolidated because of an in- 
crease in the effective length of flow path* Probably the host important 
contribution of this study is the realisation of the Importance of the 
Koeeny-Carsum equation end its modifications when applied to problems of 
fluid flew through porous media* It holds great possibilities, when sub- 
ject to the proper modification. In Its application to studies of this 
type and gives promise of eventual solution of the m ary problems of fluid 
flow still facing the petroleum industry* 
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TABLE I 

F* Lika’s 1934 at » (K is fxpari laantsl) 
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